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Abstract 

We study the geodesic motion of massive and massless test particles in the background of a 
particular class of multiple charge black holes in gauged supergravity theories in D = 4. We have 
analysed the horizon structure along with the nature of the effective potentials for the case of four 
equal charges. In view of the corresponding effective potentials, we have discussed all the possible 
orbits in detail for different values of energy and angular momentum of the incoming test particles. 
The periods for one complete revolution of circular orbits and the advance of perihelion of the 
planetary orbit have also been investigated in greater detail for massive test particles. We have 
also discussed the time period of unstable circular motion and cone of avoidance of massless test 
particles in detail. All the corresponding results obtained for massive and massless test particles 
are then compared accordingly. 
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I. INTRODUCTION: 


Over the last century, the Einstein’s theory of gravitation i.e. General Relativity (GR) 
has been extremely successful to understand various observational facts like gravitational 
redshift, the precession of Mercury’s orbit, the bending of light etc. On the other hand, 
various black hole spacetimes with or without rotation were obtained as exact solutions of 


Einstein held equations in GR namely Schwarzschild metric . Other black hole spacetimes 
including charges and rotations have also been discovered. Though GR has enjoyed the great 
success, but it is still not a complete theory to understand the physics at sufficiently smaller 
ength scale e.g. near the spacetime singularity which arises in case of gravitational collapse 
In the vicinity of the spacetime singularity, the quantum effects should be taken care 
of .enously and .ecenfly afnng fheoay haa become the paonaaing caadijfpe foa fbe aaaae 


purpose which comprises gravity in the frame work of quantum theory m- 
In continuation pursuit of the search for a quantum theory of gravity, the gauged supergrav¬ 
ity theories js ill have also captured considerable attention in recent times. In such models 
of gravity, the maximally supersymmetric gauged supergravity is realized as truncation of 
string theory or U = 11 dimensional M theory compactihed on a sphere where the gauge 
group is the isometry group of the sphere. In particular, in AdS/CFT correspondance 12| . 
this gauge group becomes the R-symmetry group of the boundary GET which is related to 
the string theory on anti-de Sitter space. Black hole solutions in the AdS sector have drawn 
much attention to understand strongly coupled gauge theory at hnite temperature on the 
GET side . For the present work, we consider four dimensional charged black holes 

in = 8^, gauged supergravity upto four charges. The construction of these black holes 
have been explicitly performed in Ref. 


15| and according to 


16| . such charged black hole 


congurations are termed as R-charged black holes. In gauged N = 8 supergravity models, 
the bosonic part of the complete Lagrangian has a negative cosmological constant A pro¬ 
portional to the square of the gauge coupling constant^ g [l^ and the black hole solutions 


are asymptotically AdS. 
Motivated from the study of mo 
of various black holes in GR 


ion of massless and massive test particles in the background 


18l-l35| and other alternative theories of gravity like string the- 


^ TV is the number of supermultiplates. 

2 A ~ —g^, where g = j and I is the AdS length scale. 
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43|, we try to make an attempt to understand the geodesics of massive neutral as 


well as massless particles in the exterior region of such R-charged black holes. Such investi¬ 
gations may indeed useful to capture the effect of cosmological constant in addition to the 
charges on the motion of the test particles. The main objective of this paper is to study all 
the possible orbits of test particles in these class of spacetimes. We analyze the equivalent 
one dimensional effective potentials in order to study the geodesics equations of motion of 
the test particles in R-charged black holes background. On the other hand, we also focus our 
study on the various parameters (mass parameter, multiple charge parameters and gauge 
coupling constant) of the solutions and present various possible orbits by tuning those pa¬ 
rameters. Due to the presence of four non zero charges, the form of the one dimensional 
effective potentials are little complicated in our case. We mainly concentrate on the numer¬ 
ical approach without performing the analytic solutions 
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-24|. The form of potentials can 


be retraced back to the known form in the limit of zero charges which otherwise takes quite 
nontrivial form. 

This paper is organised as follows. In the next section, we discuss about the spacetime of our 
interest and the horizon structure of these particular class of R-charged black holes 15| with 
multiple charges, whereas the ungauged situation (i.e. 5 ^ = 0) is same as the Schwarzschild 
black hole with flat asymptotes in GR. In section III, we derive the geodesics equation and 
obtain the form of the effective potential. In section IV, we study geodesics equation of 
motion for massive neutral particles (i.e. the timelike geodesics) under one dimensional 
effective potential by considering different values of the parameters involved in. The radial 
and non-radial geodesics are studied analytically and otherwise numerically. Using the ef¬ 
fective potential techniques, the motion of test particles and the structure of corresponding 
orbits are discussed in greater detail. In order to have a complete analysis of the geodesic 
motion in the background of the abovementioned black hole spacetime, we also study the 
motion of massless test particles (i.e. null geodesics) in section V. Finally, we conclude and 
summarize our results. 
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II. THE SPACETIME WITH MULTIPLE CHARGES 


We consider the following spacetime metric of fonr electric charge black hole solntion in 
D = 4, N = 8 ganged snpergravity, 

ds^ = -{H,H2HsH^)-^/^fdt^ + {H,H2H3H^y/\r^dr^ + r^dQl,), (2.1) 


where 

= l + f = k-^ + 2gh\H,H2H,H,). (2.2) 

Here g is the gange conpling and /i is the non-extremal parameter like mass term in pnre 
Schwarzschild black hole, k can take three valnes i.e. 1, 0 and —1, however in the present 

stndy, we consider only the case of black hole solntions with k = 1. The l3a{oi = 1.4) in 

eqs. m and fl2.2p parametrizes the fonr electric charges. 

The scalar curvatnre R for the above metric can be given as below, 


R = 


^ 1/2 




LilL 

2^n 


(2.3) 


n' ' 8""^ 

where R = HiH 2 H^H 4 . One can easily verify that for Ha = 1 (no charge solntion), the 
scalar cnrvatnre R = —8g‘^, which is exactly same as that of an Anti-de Sitter Schwarzschild 
black hole. 

Let ns hrst discnss abont the event horizons of the metric eq. fl2.ip in presence of n eqnal 
charges. So for n eqnal charges with k = 1, eqnation fl2.2p rednces to. 


H = l + -, f = 1- - + 2g'^r‘^H^, (2.4) 

r r 

with p = yusinh^/d and n = 0,1, 2, 3,4. The lapse fnnction / vanishes at the zeros of the 
eqnation given as, 

r”“^(r —/i)-f 25(^(r-I-p)” = 0. (2.5) 


The real and positive zeroes of eq. fl2.5l) are known as the horizons for the corresponding 
spacetime. Hence if a spacetime represents black hole / shonld vanish at some positive 
valne of r. For n = 1, 2, 3, there is always a cnbic eqnation for r while for n = 4, we 
have a qnartic eqnation for r. The roots of eq. (l2.5p in a closed form are really lengthy and 
complicated, so are not discnssed over here, rather we present the graphical presentation of 
the position of the horizons by imposing certain conditions among the parameters of the 
black holes. In presence of one and two eqnal charges one will always have a horizon as 
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shown in the hg([T] 2 , b, c, d) for various values of jd and g. While in presence of three and 
four equal charges, there are some constraints on the parameters involved in the metric to 
represent a black hole with at least one horizon. In presence of three equal charges, the 
constraint on parameters for the presence of horizon is given by, 

2^Vsinh®/3 < 1. (2.6) 


The three type of plots in hg([T^) or ([T]/) clearly explain the facts that black curve does not 
satisfy the bound eq. fl2.6l) with the values of parameters chosen. So it is obvious that there 
is a possibility of naked singularity in case of three equal charges if one crosses the bounds 
given in eq. (l2.6p . In case of four equal charges present, the condition of getting at least one 
horizon is given by. 


256 


< 0 , 


(2.7) 


9sinh^ g 
2 


27 

16 


A + B 

where A = sinh® jdg'^ (sinh^ fd + l)^ /r^, 
and 5 = ^ V' (sinh' fd + 2 sinh^ Id - 
The zeros of the function /(r) are presented in the plots of hg(IT^, h). In £g(IT^), one can see 
that for sufficiently small values of the constant g the number of zeros for /(r) changes. For 
certain values of Id, one will get two zeros of the function /(r) whereas smaller than that two 
zeros merge into a single zero and larger than that real zeros will disappear completely by 
leaving with a naked singularity as explained in 4^ . The same arguments are also plausible 
when the constant g varies keeping charge parameters Id hxed, as shown in fig([T}'i). 


III. THE GEODESICS EQUATIONS AND EFFECTIVE POTENTIAL 


One can simplify metric given in eq. fl2.ip by choosing d = ^, to restrict the motion of the 
particles in the equatorial plane i.e. we restrict ourselves to study of the equatorial geodesics. 
Using the following geodesic equation 


d‘^x°‘ dx^ dx'^ 

dr'^ dr dr ’ 


(3.1) 


we obtain 


t + 


f 


2n 


tr = 0, 


m\J mP 2(7 ^ 


'■W'\ ;2 




(3.2) 

(3.3) 
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FIG. 1; The lapse function /(r) for R-charged black hole: (a) and (b) represent the single charge 
case, (c) and (d) corresponds to two charge case, (e) and (f) are for three charge. Further (g) and 
(h) represent the case corresponds to four equal charges. 
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where dot denotes the derivative with respect to proper time r. 

For a timelike geodesics = — 1), 

— + + = (3.5) 


From the constant of motion, the hrst integrals of eqs. (13.2^ and (I3.4p resnlt as, 

i = EH^f-\ (3.6) 

0 = (3.7) 

where E and L are the integration constants which represent the conserved energy and 
angular momentum of the test particles respectively. 

Using eqs. fl3.6p and fl3.7p . one can now rewrite the timelike constraint eq. fl3.5p in the following 
form. 


1-0 1/FV , / 


+ - 


2\rm ^ 


E^ 

~Y' 


The effective potential for timelike geodesics may be identified from eq. fl3.8p as, 

/ / 


Veff = 


+ 1 


(3.8) 


(3.9) 


2H2 

One can reproduce the effective potential for the schwarzschild black hole with AdS asymp¬ 
tote as a limiting case when Ha = 1. It seems that under scaling of the metric components 
like / f'H~^ and —)■ one can exactly write down all the corresponding equations 

above for the case of a schwarzschild black hole in flat space However, the effective 
potential for null geodesics is given as. 


Veff = 


2rm' 


(3.10) 


The difference in the nature of effective potentials for both the cases as given in eq. (13.Op and 
eq. (l3.inp manifests itself in the structure of orbits as presented in the later sections. 


IV. NATURE OF EFFECTIVE POTENTIAL AND CLASSIFICATION OF OR¬ 
BITS FOR TIMELIKE GEODESICS 

From effective potential given in eq. (l3.9p for particles moving along the timelike geodesics, 
the equation of motion fl3.8p is investigated for two cases: radial and non-radial geodesics in 
the following subsections, respectively. 
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A. Effective potential for radial geodesics (L = 0) 


For radial geodesic (with zero angular momentum L), the effective potential given in eq. fl3.9p 
reduces to the following form, 


K// = 


/ 


2-H2 

The behaviour of effective potential is shown in fig.( 
considering all values of n ranging from 0 to 4. 


(4.1) 

for a particular set of parameters, 



FIG. 2: Variation of effective potential with r for radial geodesics with different values of n. 


If the particle is released from rest at a distance r = r*, the initial energy of the particle is 
given as. 


(r,) = 


l-^ + 2gWnir,) 


a/'^ (g) 

So from the radial equation of motion, one obtains 

1 - + 2g^r^n (r) 


in) 


which can be integrated as. 


Vn 


dr 


(4.2) 


(4.3) 


r r — 


'Vi 


E^ (n) - 


y/Hir) 


(4.4) 


In order to know the proper time experienced by a particle falling from to coordinate 
radius r, one needs to evaluate the integral given in eq. (l4.4p . Since the result in a closedform 
is not obvious here, we look for numerical solutions. The proper time given by eq. (l4.4p is 





















FIG. 3: Proper time r as a function of r with // = 1, /3 = 1, g' = 0.02 and r* = 5. Here solid, 
dotted, dashed, dot dashed and long dashed lines represent n = 0,1,2,3, 4 cases respectively and 

shows the position of event horizon. 

visually presented in fig.© which indicates that the particle falls towards the horizon in a 
finite proper time. It can also be observed from this figure that the proper time to reach 
horizon decreases with the increasing number of charges. 

B. Effective potential for non-radial geodesics (L / 0) 

Here we study the timelike geodesics for incoming test particles with nonzero angular mo¬ 
mentum. Let us try to understand orbits of the particles in the background of the R-charged 
black holes by using the effective potential 14// given by eq. fl3.9lh For r ^ (the horizon ra¬ 
dius), Veff = 0 and for r oo, 14// —t | + ^ (f “ l) sinh^(4) +sinh^(4)) ^+g‘^r‘^, 
whereas for r —>■ oo, 14// —)■ | (when g = 0). These are the asymptotic behavior of the 
potential. In fig.)© and fig.©, the 14// are given in several plots varying the charge param¬ 
eters 4 and the angular momentum L respectively keeping all the other parameters fixed. 
One can observe from figs.(jlK)-(jl^), that for a particular value of n the height of the poten¬ 
tial decreases with the increasing values of the charges. One can also notice that height will 
become small in fig.© when the values of the angular momentum are small for a particular 
value of n. 
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(a) (b) (c) 



(d) 


(e) 




FIG. 4: Effective potential for non radial geodesic with ^ = 1, g = 0.02, L = 6, represents the 
position of event horizon. Fig.(a)-(e) represent the cases of n = 0,1,2,3,4 respectively, where for 
solid line /3 = 0.1, for dashed line (3 = 1 and for dot-dashed line (3 = 1.5 (while n = 0 case does not 
contain any (3 term). 


(a) 


(b) 


(c) 
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(d) (e) 




FIG. 5: Effective potential for non radial geodesic with = 1, g = 0.02, j3 = 1, represents the 
position of event horizon. Figs.(a)-(e) represent the cases of n = 0,1, 2,3,4 respectively, where for 
ssolid line L = 2, for dashed line L = 4 and for dot-dashed line L = 6. 



FIG. 6: Effective potential for unit mass black hole (i.e. // = 1) in the presence of four equal 
charges with L = 10, /3=1, g = 0.02. 

The following orbits are allowed depending on the values of the constant E (i.e. energy of 
the incoming test particle) as plotted in hg.(l6]). 

I) E = E,-. 

(i) Here K// = E^ and r = 0 leading to a stable circular orbit at point C of hg(l6]). 

(ii) For the other possible orbit at this energy value, the test particle starts from the point 
B and falls into the singularity. Hence it is a terminating bound orbit. 
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II) E = El. 

(i) A bounded planetary orbit between points A and P (as shown in fig (l^ ). 

(ii) The other possible orbit corresponds to point D which falls into the singularity (i.e. a 
terminating bound orbit). 

III) E = E 2 . 

The possible orbits are unstable circular orbit at point E. The particle starts from point E 
and then it can go either to point G or to the singularity after crossing the horizon. 

IV) E = E 3 : 

There exists terminating bound orbit for particle crossing point J. Hence there are no fly-by 
orbits possible. 

C. Analysis of the geodesics: 

Orbit equation can be obtained by using eq. (l3.7p and eq. (l3.8p as, 

fdr\‘^ 1 , . , . 

U) 

Pn{r) = - tP'] r*. ( 4 . 6 ) 

where, / and H are given in eq. (l2.4p . The physically acceptable regions having positive 
values of Pn{r) or equivalently E'^ > V^ff for positive and real values of r, represents the 
allowed region of motion for test particles. Hence the number of positive real zeros of Pn{r) 
uniquely determine the type of particles orbit in the background of these charged black 
holes. For n = 4, the polynomial Pn{r) reduces to, 

F 4 (r) = E'^ {r + p)*^ — -|- (r -|- p)^] — r -|- (r -|- p)"^] , (4.7) 

where the mass parameter is hxed as p = 1. Note that for any value of the parameters /9, p, 
L, F, the function P 4 (r) —>■ —cx) for r —>■ ± 00 . Also, for r = 0, P^^r) = — L?' — p^) 

indicates that Pi{r) becomes positive in the domain of the parametric space where ^ > 
L? + p"^. We restrict the values of the parameters in such a way that the previous bound 
satished throughout hereafter. 

(I) For E = Ec'. 

As discussed above that for energy value E = Ec, two type of orbits are possible, one is a 
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stable circular orbit and another one is the terminating bound orbit. In fact the number of 
zeros of the function P 4 (r) uniquely characterizes the orbit structure as mentioned in 24 1 . 
Fig. ([7]) represents these cases numerically. 


(a) 


(b) 




FIG. 7: (a) Pi{r) representing the shaded allowed region for motion of test particle with a positive 
real zero at point B and a degenerate root at point C'(where the stable circular orbit exists), (b) 
Solid line represents the orbit {terminating hound orbit) for test particle starting from a distance 
at r = in the allowed region, dotted line represents the event horizon at r = rn', with L = 10, 
^= 1 , g = 0 . 02 , = 1 . 3632 . 


Circular Orbits: 


In case of circular orbit, r = tq = constant, where tq is the distance of the circular orbit 
from the singularity and hence r = 0. In order to calculate the time periods for circular 

and fl3.7p . we have 


orbits by using eqs. 


dt = 


EH{rc) 2 




Lf{rc) 

The circular orbit condition 1/= 0 can be used to express L in terms of rc as, 

A' 


L^ = 


y 


where X = - 


'c L 


+ /i 1 + 


(4.8) 


(4.9) 


-2g^U{rc) (-f+ 2 p + 2 rc), 

3; = -H-l(rc) + ^ + PhItc) ~ Tfffe) + rcH{rc) + ip) ’ 

p = /rsinh^(/3), H = 1 + ^ and H = H'^. 

As radial velocity r vanishes for circular orbit, it provides another condition for this orbit 
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as = Veff{rc). Substituting this in eq. fl4.8p with At = % and A0 = 27r for one period 
then we obtain, 


T, = 


where Tt^sch = ^ 

eqs. (l3.7p and (14.81) as, 




V}/‘^{rc)\ 


1/2 


.‘^f{rc) \rl L‘^{rc) J\ \rc + p 


rc 


T, 


t.sch') 


(4.10) 


Similarly time period in proper time can also be obtained using 


T. = 


prc f 

.2/(rc) Vc L‘^{rc) 


prc 


^rc 


1/2 


1/2 


rc 


rc +p 


27rr^ 


’'°1.2/(ro) L2(rp)'/J C.ok- (4.11) 

Now we comment on the values computed for time periods Tt and Tr respectively from 
eqs. fl4.10p and fl4.1ip . One can also easily calculate the numerical values of time periods 
for a dehnite set of parameter and a comparison of Tt and Tj. is presented in table HI In 


n 

rc 

Tt 

Tr 

0 

3.7871 

2.4851 

41.67r 

1 

4.3042 

1.9965 

44.427r 

2 

4.8169 

1.4322 

51.207r 

3 

5.3145 

1.0887 

58.267r 

4 

5.7911 

0.8640 

65.437r 


TABLE I: Comparative view of the time periods for different number of charges for fi = 1, g = 0.02 
and /3 = 1. 


comparison, one can see the effect of charges on the time periods and also on the radius of 
circular orbits, which significantly differ from the corresponding values of the Schwarzschild 
black hole ^ in GR. 

(II) For E = Ei: 

In this case, Pii^r) has three positive real zeros. The region JJ in the fig.dHt) bounded 
between two roots of P^i^r) at and rp represents the bound orbit as has been discussed 
earlier from the effective potential. Another region I in the £g.([Hti') gives us the bound orbit 

^ for Schwarzschild black holes rc = 3, Tt = 2.5375 and E = SI.IStt, for g. = 1. 
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where particle should start form the point td, shown in hg.([8)D) and crosses the horizon at 
th and then falls in to the singularity. The trajectories of the particles have been simnlated 
in figs.dHt) and (jHli). 




(0 



(d) 



-5 0 5 10 15 20 25 


FIG. 8 : (a)&:(c) PA{r) representing the shaded allowed region for motion of test particle with three 
real positive zeros at points D,P,A, (b) Solid line represents the orbit for test particle starting 
from point D in the allowed region and dotted line represents the event horizon, (d) Asterisk curve 
represents the orbit of test particle bounded between points A and P while dotted line represents 
the event horizon, ; with L = 10, /3=1, g = 0.02, = 1.75, n = 4. ((a) & (c) represent two 

magnified parts of a single plot). 


15 












Advance of perihelion for planetary orbits: 

We use the elementary derivation of the advance of perihelion of a planetary orbit for the 
Schwarzschild solution as in 45[. The advance of the perihelion for our case is obtained 
by comparing a Keplerian ellipse in a Lorentzian coordinates with the spacetime of interest 
given in eq. fl2.ip . The relevant relation communicating the two ellipse is the areal constant 
of Keplers second law. 

The following transformation of the coordinates, r and t, in the binomial approximation can 
be obtained for the spacetime used in present study, given by eq. fl 2 .ip . 


dt' = 


dr' = 


1 + 


(2/r + np) 
4r 

(2/i + np) 
4r 


np 


+ g r [l-\ -r-— 


np 


np 


g r (1 H-) (r-— 


np 


dt, 


dr. 


(4.12) 


(4.13) 


Two elliptical orbits are considered, one the classical Kepler orbit in r, t space and other for 
given spacetime in an r', t' space. In the Lorentz space one has. 


and hence the Kepler’s second law 


dA = 


dA 


rdrdcf), 


2 ^ 4 ^ 


dt 2^ dt 


Similarly in given spacetime one has. 


dA' = / TZdr'dcj), 


(4.14) 


(4.15) 


(4.16) 


where r' is given in eq. (l4.13p and the radial function TZ = r(l + ^) 4 . The binomial approx¬ 
imation of the radial function TZ is 




4r / 


(4.17) 


Therefore, applied wherever necessary the binomial approximation and the transformation 
of coordinate, we obtain the following equation after performing the integration fl4.16p by 
using eqs. fl4.13p and fl4.17p 


dA' 

dt' 


2 r 


1 + 20/i) - 


P_ 

2 

1 /3np -|- 6/i n^p‘^g‘^ 


64 


■(9np -I- 2/x) 


2 / lOnp -|- 5/i Sdn^p^g'^ 


P9 ( 


64 


+ ... 


d(j) 
- dt 


(4.18) 
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For a single orbit, 


A0' = 


^A(/)=27r 


1 + 20/i) - 

lo d4 


l/3np + 6/i nV^%n , o ^ 

-e^(9"p + 2^) 


P9 


2 / lOnp + 5/i 35n^p^g‘^ 


64 


+ .. 


d(f). 


(4.19) 


The polar form of an ellipse is given by 

P = 


I 


(4.20) 


1 + e cos (j) 

where e is the eccentricity and I is the semi-latus rectum. Therefore applying the binomial 
approximation alongwith the integration fl4.19p and using the eq. fl4.20p . we obtain the 
perihelion shift as 


671(1 Simp 571 (i 




(np + i)/- 


9 2 

Tin p 
8 




s' + 


(4.21) 


The results we have derived for the perihelion shift of the massive test particle in the space- 
time of a R-charged black holes where all n (where n = 0,1, ...,4) charges are equal. One 
can see that the corrections have appeared in eq. (l4.2ip due to the charges and the gauge 
coupling constant g apart form the term depends on mass parameter p. It is also straight 
forward to check that the form of the perihelion shift equation fl4.21l) nicely goes back to the 
results obtained for Schwarzschild AdS (i.e. n = 0, here) presented in 46l|^. 

(Ill) For E = E2: 


At this energy value of the incoming particles, the two zeros of the polynomial P 4 ,{r) in 
the previous plots (hg.dH])) merge where P^i^r) = 0 and also -^P^{r) = 0 at r = rp. The 
representative plots have been shown in £gs.(l9k) and (Eb). The radius rp is the radius of 
the circular orbit which is unstable. In £g.(j9b), a polar plot of the timelike geodesic is given 
for particle arriving from a hnite radius rp and plunges in to the black hole singularity. 
Whereas in £g.([9li) another possible orbit has shown where the particles coming from tq 
and having an unstable circular orbit at r = r^. 


^ The next order terms are discarded in I as given in equation (I4.2ip . 
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60 - 



FIG. 9: (a)&:(c) P 4 (r) representing the shaded allowed region for motion of test particle with two 
real positive zeros at points F and G, (b) Solid line represents the terminating bound orbit for test 
particle starting from point F to the singularity and dotted line represents the event horizon, (d) 
Asterisk curve represents the unstable eircular orbit of test particle while dotted line represents 
the event horizon; with L = 10, /3=1, g = 0.02, = 2.0366, n = 4. ((a) & (c) represent two 

magnified parts of a single plot). 

(IV) For E = Eg: 

At this energy value E = E^, the degenerate root ai r = rp will disappear leaving only 
a single root at r = rj as shown in fig. ffTOb ). Again we have numerically simulated the 
path of the particles carrying energy E = E^ is shown in hg. flTOh ) represents a terminating 
bound orbit. 
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(a) 


(b) 



(a) 



FIG. 10: (a) & (b) Pi{r) representing the shaded allowed region for motion of test particle with 
one real positive zeros at point J, (c) Asterisk curve represents the terminating bound orbit of test 
particle while dotted line represents the event horizon; with L = 10, /?=!, g = 0.02, E'^ = 2.5, 
n = 4. ((a) & (b) represent two magnified parts of a single plot). 

In order to have a detailed view of the geodesic motion in the background of R-charged 
Black Holes, the Null geodesics are discussed in the next section. 


V. NATURE OF EFFECTIVE POTENTIAL AND CLASSIFICATION OF OR¬ 
BITS FOR NULL GEODESICS 

A. Radial Geodesics 

Since the radial geodesics are the trajectories followed by zero angular momentum test 
particles(i.e. L=0) in given spacetime geometry, the effective potential given by eq. fl3.10p 
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for radial null geodesics vanishes i.e. 14// = 0 and consequently the radial equation of 
motion reduces to, 


<ir 


(6.1) 


where positive sign corresponds to the outgoing test particles. The integration of eq. (l5.ip 
leads to, 


r = ±Et + To, 


(5.2) 


where tq is the integration constant corresponding to the initial position of the test particle. 
Eq. fl5.2l) shows that in terms of proper time, the radial coordinate depends only on the 
constant energy value E. 

Using eq. (l3.6p and eq. (l5.ip . the radial equation of motion can also be obtained in terms of 
coordinate time t (for n = 4) as. 




(5.3) 


where / is given in eq. (l2.4p . On integrating eq. (l5.3p . the radial coordinate r can be obtained 
in terms of the coordinate time t. For the present study, the eq. (15.31) is solved numerically 
for an incoming test particle (i.e. consideration of the negative sign in eq. (15.31) 1. The radial 
motion of a massless test particle with coordinate time and proper time is presented in 
fig- dm)- One can observe from Eg. (ITTl) that a massless test particle crosses the event horizon 
in a hnite proper time while it approaches asymptotically to the event horizon in terms of 
coordinate time. 
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FIG. 11: Radial motion of an incoming massless test particle for proper time, r (dashed curve) 
and the coordinate time, t (solid curve) with ^ = 1, E' = 1, /3 = 1, n = 4, s' = 0.02; where rn ( 
verticle dot-dashed line) and represent the position of event horizon and the initial position of 
test particle respectively. 

B. Non-radial Geodesics 



FIG. 12: Effective potential for non-radial null geodesics with /x = 1, L = 6, (^ = 0.02, n = 4, 
where for solid line j3 = 0.1, for dashed line j3 = 0.5, for dot-dashed line (5 = 1, ru represents the 
position of event horizon. 
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In this section, we study the null geodesics for incoming test particles with non zero angular 
momentum. Variation in the effective potential with charge parameter (3 is shown in £g. flT^ 
while all other parameters are hxed. One can observe from fig. fll2p . that height of the 
effective potential decreases with increasing the value of charges for a particular value of n. 

The following orbits are allowed depending on the values of the constant E (i.e. energy of 



FIG. 13: Effective potential of a unit mass black hole (i. e. p = 1) for a massless test particle in 
the presence of four equal charges with L = 10, /3=1, g = 0.02. 

the incoming test particle) as shown in hg. flT^ . 

I) E = Eq: 

A particle starting form point M drops into the singularity, forming a terminating bound 
orbit. 

II) E = Ei: 

(i) A fly-by orbit is present for the particle coming from inhnity, which turns at point P as 
shown in fig. ([13]). 

(ii) For the other possible orbit at this energy value, the test particle starts from the point 
Q and falls into the singularity. Hence it is a terminating bound orbit. 

III) E = E 2 : 

A particle starting from inhnity follows an unstable circular orbit at point R which hnally 
drops into the singularity as depicted in hg (lT^ . Hence starting from point R, it follows a 
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terminating bound orbit. 
IV) E = Es: 


A particle starting from the infinity finally drops into the singnlarity after crossing the 
horizon. Hence it is a terminating escape orbit. 


C. Orbit Analysis for null geodesics: 


Using the constraint for nnll geodesics (i.e. the corresponding orbit equation reads 

as, 

here 

p^nuii^r) = rh (5.5) 

where, / and H are given in eq. fl2.4p . Again the physically acceptable regions having positive 
values of Pn^^^\r) or equivalently E"^ > Vgff for positive and real values of r, represents the 
allowed region of motion for massless test particles. The number of positive real zeros of 
p^nnZ/(^) ^] 2 iquely determine the type of orbits in the given background. Let us now analyse 
all the possible types of orbits for massless test particles. Hence, for n = 4 and /i = 1 the 
polynomial reduces to, 

p^ull^r^ _ p2 ^^4 _ p2 1-^2 _ ^ 2^2 ^ 


(I) For El = Eq: 

There exists an interesting type of geodesics for massless test particles with definite angular 


momentum, specihcally known as Cardioid type Geodesics From eq. fl5.4p . it follows 

that with specific energy value of E"^ = 2g‘^L‘^, there also exist cardioid type geodesics for 
R-charged black holes. For this special value of energy, the eq. fl5.4p reduces to, 

(S) ( 5 - 7 ) 


An elemental integration of eq. (l5.7p leads to. 


r{(j)) = ^{1 + cos (j)). (5.8) 

which is the equation of the cardioid. This solution is independent of the angular momentum 
of incoming test particle and depends only on the mass of black hole. 
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(II) For E = Ei: 

The numerical solution of orbit equation fl5.4p is shown in hg. flTT|) . Here one may notice 

(a) (b) 




(c) 



FIG. 14: (a) represents the shaded allowed region for motion of a massless test particle 

with two real positive zeros at point P and Q, (b) Solid line represents the orbit of a massless test 
particle in region II, (c) Solid line represents the orbit of a massless test particle in region I while 
dotted line represents the event horizon; with L = 10, /?=1, g = 0.02, = 1, n = 4. 

the presence of a fly-by orbit with turning point at P and a terminating bound orbit starting 
from point Q as depicted in hg. flT^ . 
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(Ill) For E = E 2 : 

The trajectory of a massless test particle with energy E = E 2 is presented nnmerically in 
hg. flTHll . The particle starting from inhnity forms an nnstable circular orbit at point R, 
which finally terminates into the singularity. 

W (b) 




(c) 



FIG. 15: (a) P 4 ^^^\r) represents the shaded allowed region for motion of a massless test particle 
one real positive zero at point R, (b) Solid line represents the orbit of a massless test particle in 
region II, (c) Solid line represents the orbit of a massless test particle in region I while dotted line 
represents the event horizon; with L = 10, /3=1, g = 0.02, = 1.736, n = 4. 

The circular orbit condition (i.e. = 0) for null geodesics is reduces to, 

2r^ — (2p + 3/i)r + /ip = 0, (5.9) 
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where K//(r) is given in eq. fl3.10p with n =4, p = /rsinh^/9. Solution of eq. fl5.9p describes 
the radius of unstable circular orbit and is given as, 



where p = /isinh^ and eq. (l5.9p has a real solution for (2p + 3/r)^ > 8/ip. The larger root of 


eq. fl5.9p locates the position of unstable circular orbit, which again has the minimum value 


for the condition (2p + Zpf = 8/ip as. 


1 (2p + 3p). 


(5.11) 



mm 


It is clear from eq. fl5.10p that radius of unstable circular orbit does not depend on the 


coupling constant g. In the absence of charges (i.e. p ^ 0), eq. fl5.10p and eq. fl5.1ip both 


reduces to the Schwarzschild case. Hence due to the presence of charges, the radius of 
unstable circular orbits is increased. 

The Time Period: The time period for unstable circular orbits can be calculated for proper 
time as well as coordinate time. The expression for time period for proper time can be 
obtained by integrating eq. fl3.7p (with 0=27r for one time period) as. 



(5.12) 


Time period in coordinate time can be obtained by combining eq. fl3.6p and eq. fl3.7p and 
integrating afterwards as follows. 



(5.13) 
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re 


Tt 

Tr 

0 

1.5000 

5.182271 

6.236671 

1 

1.6477 

7.090471 

7.361771 

2 

1.8694 

9.447271 

12.15371 

3 

2.1906 

12.18871 

19.98271 

4 

2.6174 

15.17971 

31.97971 


TABLE II: Comparison of the time periods for unstable circular orbits of massless test particles 
with different number of charges (i.e. re) for n = 1, g = 0.02 and (3 = 1. 


Cone of Avoidance: 

The cone of avoidance can be defined at any point 
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47| and the light rays included 


in the cone must necessarily cross the horizon and get trapped. If denotes the half-angle 
of the cone directed towards the black hole at large distances, then 


here 


1 dr 


if = 


,/l/2 


(5.14) 


(6.16) 


where / is given by eq. (l2.4p . The radial function TZ described in eq. (l4.16p and eq. (l4.17p is 
given as 77 = {r + p) for re = 4. Eq. fl5.15p describes the element of proper length along the 
generator of the cone. Hence, 


cot = 


1 dr 
rf^C d(j) 


(5.16) 


Now using the eq. fl5.4p . one can obtain 

tan'^ = 

From eq. fl5.17p . it follows that. 




_E\r +pY - L^frf 


TT 


1/2 


(5.17) 



V' = 2 for r = rc, 
'0 = 0 for r = ru, 


1 -|- 6p^ -|- 2(/2(4rp -|- r^) 


for r » 1, 


(5.18) 


— 2g'^) (6p2 + 4rp -j- r^) 
where, rc denotes the radius of circular orbit and ru denotes the radius of event horizon. 
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(IV) For E = E 3 : 

The path of a massless test particle with energy E = is simulated numerically in hg. flT6|l 


W (b) 




FIG. 16: (a) P/E'^^\r) represents the shaded allowed region for motion of massless test particle no 
real positive zero, (b) Solid line represents the orbit of massless test particle in the allowed region 
while dotted line represents the event horizon; with L = 10, /3=1, g = 0.02, = 2, n = 4. 

which shows the presence of a terminating escape orbit. 


VI. SUMMARY AND CONCLUSIONS 

In this article, we have investigated geodesic motion of massive as well as massless test 
particles around a particular class of R-charged black holes in N = S, D = 4 gauged 
supergravity theory. Some of the important results are summarised below: 

• It is observed that in R-charged black hole spacetimes, the horizon is always present 
for single and two charges, while it occurs only for a particular condition in case of 
charges more than two. 

• Different types of orbits (such as terminating bound, planetary, stable and unstable 
circular orbits) are present for incoming massive test particles, corresponding to their 
initial energy. No fly-by orbits are observed in any case for massive test particles. 
Fly-by and terminating escape orbits are however present for the case of massless test 
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particles unlike the case of massive test particles. There exist no stable circular orbits 
for a massless test particle. An interesting type of orbit for a massless test particle 
namely, Cardioid type geodesic is found to depend only on the mass of the black hole. 

• Radius of stable circular orbits for massive test particles and radius of unstable circular 
orbits for massless test particles are found to increase with the enhancement in the 
number of charges. We have also computed the time periods both in coordinate 
and proper time to visualise the changes with the number of charges. The radius of 
unstable circular orbit for a massless test particle is found to be independent of the 
gauge coupling constant. 

• Advances of perihelion in planetary orbits is calculated for massive test particles by 
using the Keplerian orbit method and the corrections arising from non zero charge 
parameters apart from the gauge coupling constant are obtained. 

• Cone of avoidance for a massless test particle is calculated and it is observed that at 
large distances, it depends on both the charge and gauge coupling constant. 

All the results obtained in the present study reproduces the results manifestly cor¬ 
responding to well known spacetimes such as Schwarzschild and Schwarzschild AdS 
black hole spacetimes in the prescribed limits respectively. 
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